3 results
Multi-species turbulent mixing under supercritical-pressure conditions: modelling, direct numerical simulation and analysis revealing species spinodal decomposition
- Enrica Masi, Josette Bellan, Kenneth G. Harstad, Nora A. Okong’o
-
- Journal:
- Journal of Fluid Mechanics / Volume 721 / 25 April 2013
- Published online by Cambridge University Press:
- 19 March 2013, pp. 578-626
-
- Article
- Export citation
-
A model is developed for describing mixing of several species under high-pressure conditions. The model includes the Peng–Robinson equation of state, a full mass-diffusion matrix, a full thermal-diffusion-factor matrix necessary to incorporate the Soret and Dufour effects and both thermal conductivity and viscosity computed for the species mixture using mixing rules. Direct numerical simulations (DNSs) are conducted in a temporal mixing layer configuration. The initial mean flow is perturbed using an analytical perturbation which is consistent with the definition of vorticity and is divergence free. Simulations are performed for a set of five species relevant to hydrocarbon combustion and an ensemble of realizations is created to explore the effect of the initial Reynolds number and of the initial pressure. Each simulation reaches a transitional state having turbulent characteristics and most of the data analysis is performed on that state. A mathematical reformulation of the flux terms in the conservation equations allows the definition of effective species-specific Schmidt numbers $(\mathit{Sc})$ and of an effective Prandtl number $(\mathit{Pr})$ based on effective species-specific diffusivities and an effective thermal conductivity, respectively. Because these effective species-specific diffusivities and the effective thermal conductivity are not directly computable from the DNS solution, we develop models for both of these quantities that prove very accurate when compared with the DNS database. For two of the five species, values of the effective species-specific diffusivities are negative at some locations indicating that these species experience spinodal decomposition; we determine the necessary and sufficient condition for spinodal decomposition to occur. We also show that flows displaying spinodal decomposition have enhanced vortical characteristics and trace this aspect to the specific features of high-density-gradient magnitude regions formed in the flows. The largest values of the effective species-specific $\mathit{Sc}$ numbers can be well in excess of those known for gases but almost two orders of magnitude smaller than those of liquids at atmospheric pressure. The effective thermal conductivity also exhibits negative values at some locations and the effective $\mathit{Pr}$ displays values that can be as high as those of a liquid refrigerant. Examination of the equivalence ratio indicates that the stoichiometric region is thin and coincides with regions where the mixture effective species-specific Lewis number values are well in excess of unity. Very lean and very rich regions coexist in the vicinity of the stoichiometric region. Analysis of the dissipation indicates that it is dominated by mass diffusion, with viscous dissipation being the smallest among the three dissipation modes. The sum of the heat and species (i.e. scalar) dissipation is functionally modelled using the effective species-specific diffusivities and the effective thermal conductivity. Computations of the modelled sum employing the modelled effective species-specific diffusivities and the modelled effective thermal conductivity shows that it accurately replicates the exact equivalent dissipation.
Modelling of subgrid-scale phenomena in supercritical transitional mixing layers: an a priori study
- LAURENT C. SELLE, NORA A. OKONG'O, JOSETTE BELLAN, KENNETH G. HARSTAD
-
- Journal:
- Journal of Fluid Mechanics / Volume 593 / 25 December 2007
- Published online by Cambridge University Press:
- 23 November 2007, pp. 57-91
-
- Article
- Export citation
-
A database of transitional direct numerical simulation (DNS) realizations of a supercritical mixing layer is analysed for understanding small-scale behaviour and examining subgrid-scale (SGS) models duplicating that behaviour. Initially, the mixing layer contains a single chemical species in each of the two streams, and a perturbation promotes roll-up and a double pairing of the four spanwise vortices initially present. The database encompasses three combinations of chemical species, several perturbation wavelengths and amplitudes, and several initial Reynolds numbers specifically chosen for the sole purpose of achieving transition. The DNS equations are the Navier-Stokes, total energy and species equations coupled to a real-gas equation of state; the fluxes of species and heat include the Soret and Dufour effects. The large-eddy simulation (LES) equations are derived from the DNS ones through filtering. Compared to the DNS equations, two types of additional terms are identified in the LES equations: SGS fluxes and other terms for which either assumptions or models are necessary. The magnitude of all terms in the LES conservation equations is analysed on the DNS database, with special attention to terms that could possibly be neglected. It is shown that in contrast to atmospheric-pressure gaseous flows, there are two new terms that must be modelled: one in each of the momentum and the energy equations. These new terms can be thought to result from the filtering of the nonlinear equation of state, and are associated with regions of high density-gradient magnitude both found in DNS and observed experimentally in fully turbulent high-pressure flows. A model is derived for the momentum-equation additional term that performs well at small filter size but deteriorates as the filter size increases, highlighting the necessity of ensuring appropriate grid resolution in LES. Modelling approaches for the energy-equation additional term are proposed, all of which may be too computationally intensive in LES. Several SGS flux models are tested on an a priori basis. The Smagorinsky (SM) model has a poor correlation with the data, while the gradient (GR) and scale-similarity (SS) models have high correlations. Calibrated model coefficients for the GR and SS models yield good agreement with the SGS fluxes, although statistically, the coefficients are not valid over all realizations. The GR model is also tested for the variances entering the calculation of the new terms in the momentum and energy equations; high correlations are obtained, although the calibrated coefficients are not statistically significant over the entire database at fixed filter size. As a manifestation of the small-scale supercritical mixing peculiarities, both scalar-dissipation visualizations and the scalar-dissipation probability density functions (PDF) are examined. The PDF is shown to exhibit minor peaks, with particular significance for those at larger scalar dissipation values than the mean, thus significantly departing from the Gaussian behaviour.
Direct numerical simulations of supercritical fluid mixing layers applied to heptane–nitrogen
- RICHARD S. MILLER, KENNETH G. HARSTAD, JOSETTE BELLAN
-
- Journal:
- Journal of Fluid Mechanics / Volume 436 / 10 June 2001
- Published online by Cambridge University Press:
- 22 June 2001, pp. 1-39
-
- Article
- Export citation
-
Direct numerical simulations (DNS) are conducted of a model hydrocarbon–nitrogen mixing layer under supercritical conditions. The temporally developing mixing layer configuration is studied using heptane and nitrogen supercritical fluid streams at a pressure of 60 atm as a model system related to practical hydrocarbon-fuel/air systems. An entirely self-consistent cubic Peng–Robinson equation of state is used to describe all thermodynamic mixture variables, including the pressure, internal energy, enthalpy, heat capacity, and speed of sound along with additional terms associated with the generalized heat and mass transport vectors. The Peng–Robinson formulation is based on pure-species reference states accurate to better than 1% relative error through comparisons with highly accurate state equations over the range of variables used in this study (600 [les ] T [les ] 1100 K, 40 [les ] p [les ] 80 atm) and is augmented by an accurate curve fit to the internal energy so as not to require iterative solutions. The DNS results of two-dimensional and three-dimensional layers elucidate the unique thermodynamic and mixing features associated with supercritical conditions. Departures from the perfect gas and ideal mixture conditions are quantified by the compression factor and by the mass diffusion factor, both of which show reductions from the unity value. It is found that the qualitative aspects of the mixing layer may be different according to the specification of the thermal diffusion factors whose value is generally unknown, and the reason for this difference is identified by examining the second-order statistics: the constant Bearman–Kirkwood (BK) thermal diffusion factor excites fluctuations that the constant Irwing–Kirkwood (IK) one does not, and thus enhances overall mixing. Combined with the effect of the mass diffusion factor, constant positive large BK thermal diffusion factors retard diffusional mixing, whereas constant moderate IK factors tend to promote diffusional mixing. Constant positive BK thermal diffusion factors also tend to maintain density gradients, with resulting greater shear and vorticity. These conclusions about IK and BK thermal diffusion factors are species-pair dependent, and therefore are not necessarily universal. Increasing the temperature of the lower stream to approach that of the higher stream results in increased layer growth as measured by the momentum thickness. The three-dimensional mixing layer exhibits slow formation of turbulent small scales, and transition to turbulence does not occur even for a relatively long non-dimensional time when compared to a previous, atmospheric conditions study. The primary reason for this delay is the initial density stratification of the flow, while the formation of strong density gradient regions both in the braid and between-the-braid planes may constitute a secondary reason for the hindering of transition through damping of emerging turbulent eddies.